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Anderson localization has been studied extensively for more than half a century. However, while 
our understanding has been greatly enhanced by calculations based on a small e expansion in d = 2+e 
dimensions in the framework of non-linear sigma models, those results can not be safely extrapolated 
to d = 3. Here we calculate the leading scale-dependent correction to the frequency-dependent 
conductivity ct(uj) in dimensions d < 3 . At d — 3 we find a leading correction Recr(tj) oc |w|, 
which at low frequency is much larger than the u 2 correction deriving from the Drude law. We also 
determine the leading correction to the renormalization group /3-function in the metallic phase at 
d = 3. 

PACS numbers: 05.60 Gg, 72.15 Rn 
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Introduction: Anderson localization of quantum par- 
ticles in a random potential or of classical waves in a 
random medium has been studied intensively since the 
phenomenon was first proposed in a seminal paper by 
P. W. Anderson in 1958 Jl|. I n these systems localiza- 
tion is a consequence of interference of multiple scattering 
processes, affecting the motion of particles or waves the 
more, the stronger the disorder or the more restricted the 
available geometry is, and leading to complete localiza- 
tion of quantum particles at all energies in dimensions 
d < 2. In dimensions d > 2 a quantum phase transition 
from metallic to insulating behavior takes place as the 
disorder strength A is increased beyond a certain thresh- 
old value A c . There are convincing arguments that the 
transition is continuous and that there exists a critical 
point at A = A c , characterized by a diverging correlation 
length £ (the localization length on the insulating side 
of the transition). The behavior of physical observables 
such as the conductivity a and the dielectric function 
e, as well as £, are described in terms of power laws in 
(A — A c ), with critical exponents. It has been one of the 
goals of the theory to determine these critical exponents. 

Assuming that there is only one length scale (£) in the 
system, Abrahams et al. 0] have proposed a scaling the- 
ory for the dimensionlcss conductance g as a function of 
the length L of the sample (considering only hypercu- 
bic systems), which takes the form of a renormalization 
group equation: 
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h(g)- 



(1) 



dlnL 

The /3-function has been determined in the asymptotic 
regimes of very large (metallic) and very small (insulat- 



ing) g. On the basis of these limiting behaviors it is plau- 
sible to conclude that a metal-insulator transition occurs 
in d = 3 dimensions at g w O(l). At these intermediate 
values of g a systematic and controlled calculation of the 
/3-function is not easily possible. 

The scaling hypothesis has been justified to a large ex- 
tent by a mapping of the original problem on to an effec- 
tive field theory (nonlinear a- model) of interacting matri- 
ces, first proposed by Wegner 0]. Within that model one 
confirms the result that all states are localized in d < 2 
dimensions. In order to determine critical exponents, the 
/3-function has been calculated in dimensions d = 2 + e, 
e<<l, when the critical point is shifted to large values 
of conductance, g c ^> 1, and a loop expansion in powers 
of \ j g is feasible. This program has been carried out for 
the three main symmetry classes: orthogonal, unitary, 
and symplectic symmetry [f| . In this way the /3-function 
has been found, e.g., in the orthogonal case up to terms of 
fourth order in 1/g with coefficients calculated to linear 
order in e 0, 0] ■ 

In the unitary case the leading correction terms to the 
/3-function have been calculated [H : 
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The latter result shows that in the limit of e — > the term 
linear in 1/g is absent. As a consequence, the critical 
exponent of the conductivity turns out to be s = l/2e 
which extrapolates to s = 1/2 at d = 3, and is much 
too small compared to the value of s w 1.3 obtained in 
numerical studies Q- For a review of the scaling theory 
see Q. 
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Here we propose an alternative strategy designed to 
access the /3-function in d = 3 dimensions in the metal- 
lic regime: a direct calculation of the leading terms in 
an expansion in 1/g in perturbation theory. Here we 
present first results for the unitary case. It is well known 
that the leading (one-loop) contribution is exactly zero 
in this case as the so called Cooperon propagator ac- 
quires a mass, removing its diffusion pole. As a conse- 
quence, no scale dependent (L-dependent) contribution 
to the conductance g(L) is present in the lowest order. 
Scale dependent terms growing with L may only be gen- 
erated by infrared divergent building blocks, i.e. diffusion 
propagators. It may be shown that diagrams containing 
one diffusion pole and any other elements may be com- 
bined to cancel out exactly . The argument carries over 
to higher order terms, involving more than one diffusion 
pole (but no Cooperon poles, of course), considered in 
the diffusive regime. This property is a consequence of 
gauge symmetry. 

Early attempts to calculate the quantum correction to 
the conductivity in the unitary symmetry case in per- 
turbation theory (rather than from the nonlinear sigma 
model using dimensional regularization) were discour- 
aged by the above mentioned theorem Q. However, 
by keeping track of the contributions from larger mo- 
menta q > 1//, where / is the mean free path (ballis- 
tic regime), the correct result is recovered [9j. This cal- 
culation manifestly obeys the gauge symmetry property 
discussed above and confirms the result obtained earlier 
within the dimensional regularization scheme. We now 
generalize the method to arbitrary dimensions d, in par- 
ticular d = 3. We determine the scale-dependent contri- 
butions to the conductance in the lowest nonzero order in 
1/g and from there derive corrections to the /3-function 
in d dimensions. 

Dynamical conductivity: In order to access the scale- 
dependent terms in perturbation theory, from which the 
/3-function may be derived, it is convenient to consider 
the dynamical conductivity <t{uj) of the infinite system 
as a function of positive Matsubara frequency uj, rather 
than the conductance g{L) at uo = function of 

system length L. The conductivity of an isotropic disor- 
dered system (short range disorder) is a function of the 
electron density (characterized by the Fermi wave num- 
ber /cp), and the disorder strength (characterized by the 
mean free path I) and obeys in the scaling regime the 
scaling property, expressed in terms of the correlation 
length £ 



(3) 



where z is the dynamical critical exponent, known to be 
z = (d — 2)/d [10(. We introduce the characteristic length 
L u = (D/ui) 1 / 2 , where D{oS) = a(Lo)/ud is the diffusion 
coefficient and Vd is the density of states at the Fermi 
level. In the limit lo — » the scaling function tends to a 



constant, and therefore £ = [er(w = O)] 1 / 1 - 2- ^. The scal- 
ing law Eq. (J3j> holds in the neighborhood of the quantum 
critical point (QCP). There are subdominant corrections 
to the scaling form, the most familiar one being gener- 
ated by the Drude law (using units of electrical charge 
e = 1),' 
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with nd, m the electron density and mass, r the momen- 
tum relaxation time, and Dd = v 2 f t Jd the bare diffu- 
sion coefficient. We define the correlation length £ in the 
metallic regime as <r(0) = £ 2 ~ rf . One observes that to 
the leading order £ — r V( 2 - d ) ; anc l then lot ~ Lo£~ zd , 
in conflict with scaling. However, the dominant (scale- 
dependent) contributions to <t(lo) in the limit to — > will 
obey the scaling property, as shown below. 

A renormalization group equation for o~(lo) as a func- 
tion of lo is obtained by considering the dimensionless 
conductance <? w (£ w ) = L^~ 2 a(Lo) as a function of 



dh\g u 
d In L tJ 



(5) 



In principle, the /3-function may be obtained from a cal- 
culation of the conductance in perturbation theory in the 
disorder. The dependence on length L u may be repre- 
sented in the form of a power series in ln(L w /Lo), where 
Lq is a reference length (later to be taken equal to I). One 
thus has a double power series in the disorder parameter 
A = l/(7rfcpZ) and \o.{L u /Lq) 



9u 



and 



So (A) + 5i (A) HLJL ) + g 2 (X) ln 2 (WA,) 



flv(A) = ^2 a vn X r 



(6) 



(7) 



The /3-function may be extracted from this power series 
by taking the derivative of g u with respect to In L u and 
then putting L w = Lq : 



/%o) = -Si(A( 5o )) 
9o 



(8) 



Here we defined guj(Lo) = go, and A (go) is the inverse of 
the function go (A). 

Perturbation theory: Scale-dependent quantum correc- 
tions to the conductivity are generated by diagrams with 
infrared divergent diffusion propagators (diffusons). As 
shown by Hikami [(| the perturbation theory may be 
organized in terms of diffusons connected by "Hikami 
boxes" . At two-loop order we may distinguish diagrams 
with two and with three diffusons (see Fig. [1]), with 4- 
vertex and 6- vertex Hikami boxes /14 and hg ,respectively 
(see Fig. [2]). In addition, as shown in Ref. j9(, gauge in- 
variancc requires a family of one-diffuson diagrams h 2 to 
be added within the perturbation scheme (see Fig. \S§ . 
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FIG. 1: Diagrams contributing to the 2-loop corrections to 
the conductance. Dashed squares represent 4-vertex vector 
Hikami boxes fi4, while the hexagon represents the 6- vertex 
scalar Hikami box h^. Wavy lines denote diffusons. 



h = 




FIG. 2: Diagrams contributing to the 4-vertex vector Hikami 
box hi (top), and the 6- vertex scalar Hikami box (bottom). 
Dashed lines represent impurity scattering. 



In the lowest order, we then obtain the following con- 
tribution to cr(w): 



1 f d d q , f d d Q 
^-loopH = -J J^yDMj J^y F M>")> 

(9) 

1 
d 



(27T) 

F(q,Q,(j)=D(Q,cj) 



h 2 4 (q,Q)D(<i + Q,u J ) 



h 6 (q,Q) + h 2 (Q) 



(10) 



In this expressions we have singled out one diffuson with 





FIG. 3: The additional 2-vertex (one-diffuson) diagrams /12 
required by gauge invariance for 2-loop corrections to the con- 
ductance. 



momentum q and collected all the remaining factors in 
the integral of the function F(q, Q, lu). 

Quite generally the diffusion propagator D(q,uj) is ob- 
tained by summing the particle-hole ladder diagrams as 



D(q,u) 



2 



1-71 



imp 
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where uq is a short range impurity scattering potential, 



is the impurity density, and 



n(«, w) = y G k +q (^ + iu>)G£(E F ) (12) 

k 

Here the retarded (advanced) Green's function of elec- 
trons of mass m is given by 

cmE) = lGt(E)r = EF _ k2/ l m + t/2T , (13) 

with Ep the Fermi energy and r the single particle re- 
laxation time, defined in terms of the disorder potential 
by t^ 1 = 2-KV d n imp ul. 

We are interested in the momentum regime q < kp, 
anticipating that all momentum integrals are convergent 
and the contribution from the momenta q > kp is neg- 
ligible (kp is the Fermi wave vector). We consider the 
weak disorder regime only, where kpl 3> 1. In the limit 
of small momentum, g( < 1 the diffuson has the usual 
diffusion pole form 



D(q,u) 



1 



1 



(14) 



2-Kv d T 2 D d q 2 + uj 

In the same diffusive limit ql,Ql <C 1, one has to retain 
only the diagram (a) in the 4-vertex Hikami box and 
the diagrams (a-c) for the 6-vertex Hikami box he (see 
Fig. 2) while the single diffuson diagrams of Fig. 3 are 
not important. This yields 

hl(Q,q) = -(47r^Ai) 2 T 6 (Q + q) 2 , (15) 
h 6 (Q,q) = AKv d D dT \ (16) 
h 2 (Q) = 0. (17) 

We now substitute these results into the above expres- 
sion for the conductivity, first concentrating on the con- 
tribution from the diffusive regime 



°2-loopM 
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2-K 2 v d J (2n) 2d D d q 2 +uD d Q 
2D d {<i+Q) 2 /d 
D d (q + Q) 2 +Lj 



1 - 



(18) 



In low dimensions, d < 2, the integrals are convergent, 
yielding the scale-dependent result 



°2-loop( w ) = C d<2 v d D d 



2/{2-d) n d/(2-d) 
UV d U d 



d-2 



Cd<2£ 



2-d 



4-2d 



, (19) 
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where the constant is given by 
1 f d d x d d y 1 



Cd<2 = 



dn 2 



(2ir) 2d x 2 + ly 2 + l 



- 1 



(x + y) 



1 



(20) 



The integrals may be reduced to a single integral of a 
hypergcomctric function 



<2 



-2d-d-2 
TV 
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K l dtF { 2 - d 2 d 2^-r l 



2T(d - 3) 
' d(4 - d) 



(21) 



with Ci = -1/(24tt 2 ) and C 2 _ e = -1/(16ttV), at e < 1. 
Thus, in the limit e — > 0, we find the logarithmic scaling 
dependence 



(T. 



2-loop 
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In L„ 
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(22) 



which is in agreement with the known two-loop correc- 
tion result to the conductivity in d = 2 dimensions (here 
obtained by approaching d — 2 from below). 

The same result can be obtained directly in d = 2, 
without resorting to the e expansion. If we assume that 
q < Q, then the infrared-divergent part of Eq. ([§]) origi- 
nates from the integral over q, while the Q-integral yields 
the numeric prefactor. The diffusons and Hikami boxes 
entering the function F(q, Q,cj) should be calculated 
with an increased precision allowing for high (ballistic) 
values of Q. This requires taking into account all the 
diagrams shown in Figs. 2 and 3. Ballistic contributions 
make the integral over Q convergent leading to the same 
logarithmic correction (f22|) . 

We now turn to dimensions d > 2, for which the mo- 
mentum integrals in Eq. (|18p are not converging. In order 
to regularize these integrals we first take a logarithmic 
derivative with respect to frequency ui 



^t*oM = h + h 



where 



h = - 



D d fd d qd d Q 



1 



ir 2 is d J (2ir) 2d (D d q 2 + uj) 2 D d Q 2 + lj 

2Aj(q+Q) 2 AT 
D d (q + Q) 2 +LJ_ 



d d q d d Q oj 



dir 2 v d J {2ir) 2d D d q 2 + u 

1 £> d (q + Q) 2 



D d Q 2 +Lu[D d {q+Q) 2 +u} 2 - 



(23) 

(24) 
(25) 
(26) 
(27) 



While 7 2 converges provided d < 3, (for a further reg- 
ularization see below), I\ is still ultraviolet divergent. 
Including the contribution to the momentum integral on 
Q from the ballistic regime we have 



h 



A-K 2 u d r 2 J (2ir) d {D d q 2 + uof 



d d Q 



jF(q,Q, w). (28) 



(27T) 



We now use the identity / d d QF(0, Q, 0) = 0, which may 
be checked by direct calculation and, more generally, is 
a consequence of the gauge invariance Therefore, 
replacing F(q,Q,w) by F(q,Q,w) - F(0,Q,0) in Eq. 
(|28j) . we obtain convergent integrals in the ultraviolet, 
which allows us to use the diffusive limit expressions for 
-F(q, Q,w). Adding the contribution 7 2 , and reformulat- 
ing the integrals using the symmetry under q <H> Q, we 
then find 



dlnuj 



oopH 



'2-loop 



d In uj dir 2 v d 



d d qd d Q 



1 



(2ir) 2d D d q 2 
1 



u DdQ 2 + ui 
(1 - d/2)u 



D d (q + Q) 



D 2 q 2 Q 2 



(29) 



The scaling contribution to the conductivity in dimen- 
sions 2 < d < 3 takes the same form as in Eq. (fT9"]) where 
the constant C is now given by 



a 



d>2 



1 

dn 2 



d d x d d y 1 



1 



(2^) 2 
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x z y £ 



+ 



(x- 



1 



(30) 



The integrals may be done to give the same expression as 
for d < 2, namely Eq. ((21]) . This means that by analytic 
continuation in the complex d-plane one may pass from 
the regime d < 2 to the regime d > 2, around the singu- 
larity at d = 2. For dimensions close to integer values we 
get C 2+e = l/(167r 4 e) = -C 2 _ e , and C 3 - e = l/(967r 4 e). 

The logarithmic ultraviolet divergence in the case of 
three dimensions is cutoff by the upper limit q = q rj 
1//, leading to a logarithmic contribution 



T ti=3 
T 2-\oo 



i In. c 



967r 4 j/ 3 Li2 



48tt 



V-'J™. (31) 



Subleading corrections requiring a precise calculation of 
the contribution from the ballistic regime lead to the fac- 
tor l/£ and the constant c under the logarithm in the last 
equation. We now analytically continue from the imagi- 
nary frequency axis to the real axis by replacing the Mat- 
subara frequency uj by the real frequency O = iu. We 
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note that Rea^ op (ft) ~ |0| and Im ^_y o 3 op (ft) ~ filn 
in that case. 

It is remarkable that the correction term to Rea(fl) 
varies linearly with 57, 



Rea d=3 



loop(^) 



9-7T 



H 



(32) 



and is thus dominant at frequencies fl < f2 s , compared 
to the leading finite frequency correction ~ — cto(£2t) 2 
of the Drude law. Moreover the scale-dependent term 
is positive in contrast to the Drude correction, so that 
we predict that the conductivity has a maximum at oj = 
U s / 2. We estimate the crossover frequency as 



tt s = ^-{k F l) 5 e F . 



(33) 



At finite temperature the RG-flow is cut off by phase 
relaxation effects. This is because in the presence of in- 
teraction the number of particles at given energy is no 
longer conserved, giving rise to a phase relaxation term 

1 /t^ in the independent particle diffusion pole expression 
Eq. (|I4j) . The phase relaxation rate is temperature de- 
pendent and vanishes at T = (for calculations of 

in two dimensions see [l2, 13|). At finite T one may re- 
place £2 by f2 + i/t0 in Eq. (f31~j) . At SI <C i/r$ one then 
finds a ~ il 2 T0, which is distinguished from the Drude 
contribution by a temperature dependent prefactor grow- 
ing with decreasing temperature and is larger than the 
Drude term if < Q s . 

Renormalization group equation: We are now in a posi- 
tion to determine the leading term in a 1 /^-expansion of 
the RG-/3-function for any dimension d. For dimensions 

2 < d < 3, we then find, following the steps outlined in 
the introduction 



d In 

d In L u 



2{d~2)C d - 



(34) 



In the limit d — > 2, putting d = 2 + e, we recover the 
known result (now approaching d = 2 from above) 



1 

8tt 4 



1 



(35) 



The case d = 3 requires special attention. As dis- 
cussed above, in that case a logarithmic correction factor 
in the scale dependence is found 
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Consequently, the 
/3-function also acquires a logarithmic correction factor 



1 hi 3^ 



C'S 



48tt 4 



(36) 



The small prefactor 1/48-7T 4 makes the correction term 
small and probably insignificant in the neighborhood of 
the transition point. The latter statement anticipates 
that the constant C3 is of order unity and therefore does 



not shift the critical point much. Preliminary calcula- 
tions of the contributions to C3 appear to confirm this 
assumption. It therefore appears that our result does 
not allow a reasonable determination of the critical ex- 
ponent. It rather suggests that the /3- function has a more 
complex shape caused by higher loop contributions. 

Conclusion: In the above we presented a derivation 
of the leading scale-dependent contribution to the dy- 
namical conductivity of a disordered system in the uni- 
tary symmetry class within the model of non-interacting 
fermions. We showed how the scaling terms may be ex- 
tracted from the perturbation theory in any dimension 
d. On the basis of these results we determined the lead- 
ing term in the renormalization group /3-function in the 
regime of dimensions 2 < d < 3. It is interesting to note 
that the leading correction term in /3 at large conduc- 
tance g varies as g~ 2 , independent of dimension, with a 
logarithmic correction at d = 3. Its prefactor is rather 
small, suggesting that higher loop order terms will be 
important near the transition point (/3(<? c ) = 0). It will 
therefore not be easily possible to determine the critical 
exponent in that way. 
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